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We discuss the signature of a cosmic string wake in 21cm redshift surveys. Since 21cm surveys 
probe higher redshifts than optical large-scale structure surveys, the signatures of cosmic strings 
are more manifest in 21cm maps than they are in optical galaxy surveys. We find that, provided 
the tension of the cosmic string exceeds a critical value (which depends on both the redshift when 
the string wake is created and the redshift of observation), a cosmic string wake will generate an 
emission signal with a brightness temperature which approaches a limiting value which at a redshift 
of 2: + 1 = 30 is close to 400 mK in the limit of large string tension. The signal will have a specific 
signature in position space: the excess 21cm radiation will be confined to a wedge-shaped region 
whose tip corresponds to the position of the string, whose planar dimensions are set by the planar 
dimensions of the string wake, and whose thickness (in redshift direction) depends on the string 
tension. For wakes created at Zi + 1 = 10^ , then at a redshift of z -f 1 = 30 the critical value of the 
string tension fi is Gfj. = 6 x 10~^, and it decreases linearly with redshift (for wakes created at the 
time of equal matter and radiation, the critical value is a factor of two lower at the same redshift). 
For smaller tensions, cosmic strings lead to an observable absorption signal with the same wedge 
geometry. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Cosmic strings (see [IHll for initial work on cosmic 
strings and structure formation) cannot be the dom- 
inant source of the primordial fluctuations, however they 
can still provide a secondary source of fluctuations. Over 
the past decade, the realization has grown that many 
inflationary scenarios constructed in the context of su- 
pergravity models lead to the formation of gauge theory 
cosmic strings at the end of the inflationary phase 9]. 
Also, in a large class of brane inflation models the forma- 
tion of cosmic superstrings at the end of inflation is 
generic and in some cases (see these strings are 
stable (see also [ll| for reviews on fundamental cosmic 
strings). Cosmic superstrings are also a possible remnant 
of an early Hagedorn phase of string gas cosmology [3] . 
In all of these contexts, both a scale-invariant spectrum 
of adiabatic coherent perturbations and a sub-dominant 
contribution of cosmic strings is predicted. Hence, it is 
important to search for the existence of cosmic strings. 

In this paper we study the signature of cosmic strings 
in 21cm radiation maps (see e.g. [l5| for an in-depth re- 
view of 21cm cosmology). Observing the intensity of the 
cosmological background radiation at wavelengths corre- 
sponding to the red-shifted 21cm transition line of neu- 
tral hydrogen has several potential advantages compared 
to the currently explored windows . First of all, it probes 
the distribution of the dominant form of baryonic matter 
and is thus not sensitive to our incomplete understand- 
ing of star formation and non-linear evolution, which is 
a problem when interpreting the results of optical red- 
shift surveys. Secondly, it probes the universe at higher 
redshifts and allows us to explore the "dark ages" (the 
epoch before star formation and non-linear clustering set 



in). Related to this, it explores the distribution of mat- 
ter in a regime when the amplitude of the fluctuations 
is smaller and linear theory is a better approximation. 
Finally, 21cm surveys provide three-dimensional maps, 
a significant potential advantage over CMB anisotropy 
maps. 

Cosmic strings are known to give rise to distinctive sig- 
natures in CMB temperature anisotropy maps fld\ , CMB 
polarization maps [17,] and large-scale structure (LSS) 
maps pi^| - [23j . These distinctive signatures come from 
moving long (compared to the Hubble radius) strings 
(see Section 2). In the CMB temperature maps, these 
strings lead to line discontinuities, in the polarization 
maps to (roughly) rectangular regions with extra polar- 
ization, and in LSS maps to thin planar regions of en- 
hanced density. These signals are manifest in position 
space maps, but they become obscured when calculating 
power spectra. Hence, the lesson is to study the maps in 
position versus momentum space. 

Here, we compute the signature of a single cosmic 
string wake in 21cm emission. We find that strings 
with tensions fi somewhat below the current limits of 
Gfi = 3 X 10^^ could be detected in 21cm maps where 
they would appear as wedges in 21cm maps with either 
extra emission or extra absorption, depending on the ten- 
sion and on the specific redshift. The planar dimensions 
are set by the direction of the string and its velocity vec- 
tor, and the width of the wedges is proportional to the 
string tension. This signal must be searched for in posi- 
tion space maps. In Fourier space the distinctive phase 
information would be washed out. 

The outline of this paper is as follows: In Section 
2 we discuss how string wakes are generated and how 
these lead to distinctive signals for observations. We re- 
view gravitational accretion of matter onto cosmic string 
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wakes and compute the temperature of the HI gas in- 
side the wake. In Section 3 we then study the 21cm 
emission signal from a single cosmic string wake. The 
generalization to the case of a network of string wakes 
will be addressed in a future publication. We compute 
the brightness temperature and describe the geometrical 
structure of the signal, a structure very characteristic for 
cosmic strings. In the final section we summarize our 
results and put them in the context of other work on 
the possible detection of cosmic strings. When comput- 
ing the magnitude of the 21cm signal, we use the same 
WMAP concordance values [l^l for the cosmological pa- 
rameters as used in the review llSll. 



II. COSMIC STRINGS AND LARGE-SCALE 
STRUCTURE 

In many field theory models, the formation of a net- 
work of cosmic strings is an inevitable consequence of 
a phase transition in the early universe (for reviews on 
cosmic strings see e.g. |25l - [28| ). This network of cos- 
mic strings approaches a "scaling solution" which means 
that the statistical properties of the string network are 
the same at all times if distances are scaled to the Hub- 
ble radius [2^ . The detailed form of the scaling solution 
must be obtained from numerical simulations [20| of the 
evolution of a network of cosmic strings (see also [sO] 
for some recent analytical work). There are two compo- 
nents of the string network - firstly a network of "long" 
strings with a mean curvature radius = cit, where ci is 
a constant of order unity, and secondly a distribution of 
string loops with radii smaller than the horizon which re- 
sults from the "cutting up" of the long string network as 
a consequence of string intersections. According to more 
recent cosmic string evolution simulations, the long string 
component is more important for cosmological structure 
formation. 

Numerical simulations by different groups have clearly 
verified the scaling solution for the long string network. 
There is, however, still a large uncertainty concerning the 
distribution of string loops (the only agreement seems 
to be that the loops are less important than the long 
strings for cosmological structure formation). Hence, in 
order to obtain constraints in models with cosmic strings 
which are robust against the uncertainties in the numer- 
ical simulations, it is important to focus on signatures 
of long strings as opposed to signatures of string loops. 
The tightest current constraints on cosmic strings come 
from the shape of the angular power spectrum of CMB 
anisotropies, yielding a constraint of [31j] (see also }32i] ) 

G^i < 3 X 10-^ (1) 

where ^ is the string tension and G is Newton's constant. 
However, it is important to keep in mind that this limit is 
based on an analytical description of the scaling solution 
which contains a number of parameters which can only 



be determined from comparisons with numerical simu- 
lations and whose values thus have a substantial uncer- 
tainty. Work based on a field theory simulation of strings 
[11] gives a limit of twice the above value (see also 
for a very recent analysis of the different bounds). Note 
that limits on the cosmic string tension which come from 
gravitational radiation from string loops [SSj are sensi- 
tive both to the large uncertainties in the distribution of 
string loops, and also to back-reaction effects on cosmic 
string loops (see e.g. [1^), and are hence not robust. 

Direct limits on the cosmic string tension fi can be 
obtained by looking for specific signatures of individual 
long strings. These are limits which are insensitive to the 
parameters in the cosmic string scaling solution. Limits 
obtained from searching for the Kaiser- Stebbins signa- 
ture of a l ong string in CMB temperature maps were 
derived in [11 Hi] using WMAP data. The limits ob- 
tained were weaker than the ones from ([1]). However, the 
angular scale of the WMAP experiment is too large to 
be able to effectively search for sharp features in position 
space maps such as those predicted by cosmic strings. It 
was pointed out [39| that the string signatures for values 
of Gfi somewhat smaller than the limiting value of ^ 
should be clearly visible in smaller angular scale CMB 
anisotropy maps such as those provided by the ACT [43| 
and SPT [4l| experiments. In recent work (42| - |4^ it was 
shown that limits up to an order of magnitude tighter 
than ([Ij might be achievable using SPT data. In the fol- 
lowing, we will discuss direct signals of wakes created by 
long strings in 21cm surveys. 

As first pointed out in [ISj] and then further discussed 
in [l9l - [2^ , long strings moving perpendicular to the tan- 
gent vector along the string give rise to "wakes" behind 
the string, i.e. in the plane spanned by the tangent vec- 
tor to the string and the velocity vector. The wake arises 
as a consequence of the geometry of space behind a long 
straight string [H, HI] - space perpendicular to the string 
is conical with a deficit angle given by 

a = SirGfi. (2) 

From the point of view of an observer behind the string 
(relative to the string velocity vector), it appears that 
matter streaming by the string (from the point of view of 
the observer traveling with the string) obtains a velocity 
kick of magnitude 

Sv = AttGuVs^s ~ 4km/s(G^)6Ws7s (3) 

towards the plane behind the string. In the above, Vs is 
the velocity of the string (in units of the speed of light), 
7s is the corresponding relativistic gamma factor, and 
(G/i)6 is the value of G/i in units of 10^^. This leads to 
a wedge-shaped region behind the string with twice the 
background density (see Figure 1). 

Since the strings are relativistic, they generally move 
with a velocity of the order of the speed of light. There 
will be frequent intersections of strings. The long strings 
will chop off loops, and this leads to the conclusion that 
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FIG. 1: Sketch of the mechanism by which a wake behind a 
moving string is generated. Consider a string perpendicular 
to the plane of the graph moving straight downward. From 
the point of view of the frame in which the string is at rest, 
matter is moving upwards, as indicated with the arrows in 
the left panel. From the point of view of an observer sitting 
behind the string (relative to the string motion) matter flow- 
ing past the string receives a velocity kick towards the plane 
determined by the direction of the string and the velocity vec- 
tor (right panel). This velocity kick towards the plane leads 
to a wedge-shaped region behind the string with twice the 
background density (the shaded region in the right panel). 



the string distribution will be statistically independent 
on time scales larger than the Hubble radius. Hence, 
to model the effects of strings we will make use of a 
toy model introduced in pT] and used in most analyti- 
cal work on cosmic strings and structure formation since 
then: we divide the time interval between the time of 
equal matter and radiation and the current time to into 
Hubble expansion time steps. In each time step, we lay 
down a distribution of straight string segments moving 
in randomly chosen directions with velocities chosen at 
random between and 1 (in units of the speed of hght). 
The centers and directions of these string segments are 
also chosen randomly, and the string density corresponds 
to N strings per Hubble volume, where N is an integer 
which is of the order 1 according to the scaling of the 
string network. The distribution of string segments is 
uncorrelated at different Hubble times. Each string seg- 
ment will generate a wake, and it is the signal of one of 
these wakes which we will study in the following. 

A string segment laid down at time ti will generate a 
wake whose dimensions at that time are the following: 



CiU X tiVsJs X iTTG^tiVs^s 



(4) 



where ci is a constant of order one. In the above, the 
first dimension is the length in direction of the string, the 
second is the depth. The fact that the string wake has a 
finite depth is due to the causality constraints on density 
fiuctuations produced during a phase transition j4^] : The 
information about the formation of a string (and hence 
the information about the existence of a deficit angle) 
cannot have propagated farther than the horizon from 
the point of nucleation of the defect. It was shown in 



[49( 1 that the deficit angle quite rapidly tends to zero as 
the horizon is approached. 

The overdense region in the wake will lead to the grav- 
itational accretion of matter above and below the wake 
towards the center of the wake. In this way, the wake 
will grow in thickness. The accretion of matter onto a 
cosmic string wake was studied in 12 ll. ISOl in the case of 
the dark matter being cold, and in |50l . l5ll | in the case of 
the dark matter being hot. We are interested in the case 
of cold dark matter. As an aside, we mention that if the 
strings have lots of small-scale structure then they will 
have an effective tension which is less than the effective 
energy density [l^l • This will lead to a local gravitational 
attraction of matter towards the string, a smaller trans- 
verse velocity, and hence to string filaments instead of 
wakes. The gravitational accretion onto string filaments 
was studied in [53]. 

We will now review the computation of the width of 
the wake. We are interested in the distribution of baryons 
in the vicinity of the string. However, for t > tree the 
baryons and cold dark matter feel the same gravitational 
attraction and will thus behave in the same way - mod- 
ulo thermal velocity effects which will be discussed later - 
until the baryons undergo a shock. Thus, for the moment 
we focus on the onset of clustering of the cold dark mat- 
ter. We are interested in wakes created after the time teq 
of equal matter and radiation (there is no gravitational 
clustering of cold dark matter before that time). For 
times between teq and recombination (tree), the baryons 
are coupled to the radiation. However, for t > tree the 
baryons will rapidly fall into the potential wells created 
by the cold dark matter and thus, once again, it is le- 
gitimate to focus attention on the clustering of the cold 
dark matter. Note that wakes produced at the earliest 
time are the most numerous. 

We consider a mass shell located at a comoving dis- 
tance q above the wake. Its physical distance above the 
wake at time t is 



w{q,t) = a{t){q~tp) 



(5) 



where tp is the co-moving perturbation induced by gravi- 
tational accretion. If the wake is laid down at the initial 
time ti (with corresponding redshift Zj), then the initial 
conditions for the cold dark matter fluctuation are 



^it^) = i>{ti) = 0. 



(6) 



As a consequence of the initial wake planar overdensity 
(j{ti), the CO- moving displacement if) will begin to in- 
crease. The clustering dynamics can be studied mak- 
ing use of the Zel'dovich approximation [sj] in which the 
gravitational force is treated in the Newtonian limit. As 
reviewed recently in [l?] , we obtain 

,,,, ISttG /, n /ii \2/3 / t n2/3,2 /»n 

m = ^^^(*')(^) i-) ^0, (7) 

where to is the present time and a{ti) is the initial wake 
planar density excess. 
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The mass shell with initial comoving distance q above 
the wake "turns around" when w{q,t) = 0. At time t, 
this occurs for a value q ~ q^i (t) given by 

qmit) = ^Gfiv,j,{z, + l)-i/2i„(l)2/3 
It is easy to check that at the point of turnaround 

'^{qnut) = ^g, (9) 

and that hence the density inside the turnaround surface 
is twice the background density. 

Once a matter shell reaches its maximal distance Wmax 
from the center of the wake, it will start to collapse 
onto the wake. The infall of matter will halt as the 
shell hits other streams of matter. This will lead to a 
shock. In analogy with what can be shown analytically 
in the context of the spherical collapse model (see e.g. 
^5|) we assume that the shock occurs at approximately 
half the maximal distance. The hydrodynamical simula- 
tions of [56| confirm the applicability of this assumption. 
Since the distance at turnaround is half the width the 
shell would have without gravitational accretion onto the 
wake, the shock occurs at a distance 1/4 of that which 
the matter shell would have under unperturbed Hubble 
expansion. Hence, the average density inside the wake is 
four times the background density, a result which will be 
used several times in the computations of the following 
section. 

The evolution of the mass shell between turnaround 
and shock can be followed using (O and ([7]). The shock 
will occur when w{q,t) — [l/'i)wmax{q)- A straightfor- 
ward computation shows that at this point the velocity 
is given by 

5 z + 1 

It is this velocity which then determines the temperature 
inside the shocked region. 

The shocks will lead to thermalization of the wake. 
The gas temperature will be given by 

hBT=]^mv\ (11) 

where m is the mass of a HI atom and v is the velocity 
of the in-falling particles when they hit the shock, which 
is given by 

v = w{q{t),t), (12) 

where q{t) is the comoving distance of the shell which 
starts to collapse at the time t. 

Inserting the result ([TU]) into ((T^ and then into PT|) 
we obtain the following result for the temperature of 
the HI atoms inside the wake 

^ [20K]{Gfi)l{vsjsr^, (13) 



where in the second line we have written the result in 
degrees K and expressed G^ in units of 10~^ (and kept 
only one significant figure in reporting the final number) . 
The wake temperature is largest for wakes produced at 
the earliest times since the initial wake density is then 
the highest, and increases as time increases because more 
matter has time to accrete. 

Assuming values of Gfi = 3 x 10~^ (the current upper 
bound on the string tension), Zi — 10"^, (close to the time 
of recombination) , z+1 — 30 and (ws7s)^ = 1/3, Eq. (fT3|) 
yields Tk ~ 20K. This temperature is smaller than the 
CMB temperature ~ 82K at this redshift. As we will 
see in the following section, this leads to an absorption 
signal in the 21cm radiation. 

The formation of a cosmic string wake and the ther- 
malization which takes place inside the shocked region 
has been studied in [s^ using an Eulerian hydro code 
(S?! optimized to resolve shocks. The numerical simula- 
tions of [s^ show that the density and temperature inside 
the shocked region are indeed roughly uniform and that 
the temperature obtained agrees with the values obtained 
here using analytical approximations. 

Looking for signals from cosmic strings in 21cm surveys 
is potentially more powerful than looking in large-scale 
optical redshift surveys. This is because firstly the non- 
Gaussian signatures from strings are more pronounced at 
higher redshifts, secondly because we are directly looking 
at the distribution of the baryons, and not just the dis- 
tribution of stars, the latter being affected by non-linear 
and gas dynamics, and thirdly because the distribution 
of matter is more linear at higher redshifts and hence eas- 
ier to follow analytically. Compared to CMB and CMB 
polarization maps, 21cm surveys have the advantage of 
yielding three- rather than just two-dimensional maps, 
maps thus containing much more information. To our 
knowledge, there has been little previous work on 21cm 
emission from strings. Two exceptions are (ssj in which 
the angular power spectrum of 21cm emission from a net- 
work of cosmic strings was considered, and [soj . a recent 
study in which the correlation between 21cm emission 
and CMB signals from cosmic strings was studied. In 
contrast to these works, we are looking for direct string 
signals in position space 21cm maps. 

In the following section we will briefly summarize some 
key features of 21cm cosmology and apply the equations 
for the 21cm brightness temperature to the case of emis- 
sion from a string wake. 



III. COSMIC STRINGS AND 21CM MAPS 

Neutral hydrogen is the dominant form of baryonic 
matter in the "dark ages", i.e. before star formation. 
Neutral hydogen has a 21cm hyperfine transition line 
which is excited if the hydrogen gas is at finite temper- 
ature. Hence, we expect redshifted 21cm radiation to 
reach us from all directions in the sky, and the inten- 
sity of this radiation can tell us about the distribution of 
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neutral hydrogen in the universe, as a function of both 
angular coordinates in the sky and redshift. Thus, in 
contrast to the CMB, 21cm surveys can probe the three- 
dimensional distribution of matter in the universe (see 
for pioneering papers on the cosmology of the 21cm 
line and .15. ] for an in-depth review). 

Let us now consider the equation of radiative transfer 
along the line of sight through a hydrogen gas cloud of 
uniform temperature- in the case of interest to us this 
gas cloud is the cosmic string wake. The brightness tem- 
perature Tb(i/) at an observed frequency v due to 21cm 
emission is given by 



(14) 



where Tg is the spin temperature, T-y is the microwave 
radiation temperature, and is the optical depth ob- 
tained by integrating the absorption coefficient along the 
light ray through the gas cloud. The frequency v is the 
blue-shifted frequency at the position of the cloud cor- 
responding to the observed frequency Vq. The term pro- 
portional to Ts is due to spontaneous emission, while the 
term proportional to is due to absorption and stimu- 
lated emission. 

As explained in [15], what is of observational interest is 
the comparison of the temperature coming from the hy- 
drogen cloud with the "clear view" of the 21 cm radiation 
from the CMB. 



l + z 



Ts-T^{u) 
l + z 



(15) 



Note that the "clear view" of the CMB is hypothetical 
since even without a string wake's gas cloud the inter- 
galactic medium is partly a less dense hydrogen gas. In 
the second part of the equation above we have expanded 
the exponential factor to linear order in the optical depth. 

The spin temperature Ts is defined as the relative 
number density of atoms in the hyperfine energy states 
through ni/no — 3exp(— T^/Ts) . Here ni and no are 
the number densities of atoms in the two hyperfine states, 
and Ti, = i?io/fcs = 0.068K is the temperature corre- 
sponding to the energy splitting Eiq between these states. 

The spin temperature is determined solely by the tem- 
perature Tk of the gas in the wake, as long as UV scat- 
tering is negligible (which is true in our case). The re- 
lationship between spin and kinetic gas temperatures is 
expressed through the collision coefficients Xc which de- 
scribe the rate of scattering among hydrogen atoms and 
electrons: 



Ts 



1 



1 + XcT^/Tk 



(16) 



We will shortly discuss the numerical values of the Xc for 
a particular case of interest. 

Combining (IT5|) and (IT51) we find for the difference 
STh{i') in brightness temperature induced by the inter- 
action of the radiation with neutral hydrogen is 



(17) 



where the last factor represents the red-shifting of the 
temperature between the time of emission and the 
present time. It is important to keep in mind that the 
brightness temperature from a region of space without 
density perturbations is nonzero. It is negative since the 
temperature of gas after redshift 200 is lower than that of 
the CMB photons, the former red-shifting as {z{t) + 1)^, 
the latter as z(i) -t- 1. The optical depth of a cloud of 
hydrogen is 



3c^AiQ , hvxQ s Nhi 



4i/2 



'ksTs' 



(18) 



where Nfjj is the column density of HI. 

Up to this point, the analysis has been general. Let us 
now specialize to the case where the gas cloud is the gas 
inside the cosmic string wake. In this case, the column 
density is the hydrogen number density n^j^^ times the 
length that the light ray traversed in the cloud. This 
length will depend on the width w of the wake and the 
angle 6 that the light ray makes relative to the vertical 
to the wake so that 



Nhi = 



COS 6 



(19) 



The factor of 2 results since w is the width of the wake 
from the center. 

The line profile 4>{v) is due to broadening of the emis- 
sion line. This broadening reflects the fact that not all 
photons resulting from the hyperfine transition will leave 
the gas cloud at the same frequency. Frequency differ- 
ences are in general due to thermal motion, bulk motion 
and pressure effects. Since the pressures we are con- 
sidering are small by astrophysical standards, pressure 
broadening is negligible. Since the gas temperature inside 
the wake is not much larger than the CMB temperature, 
thermal broadening will not be important, either. This 
leaves us with the effects of bulk motion, more specifi- 
cally the expansion of the wake in planar directions. The 
line profile is normalized such that the integral of it over 
frequency is unity. 

The origin of the line broadening due to the expansion 
of the wake is illustrated in Figure 2. Let us consider a 
point on the wake for which the photons travel to us at 
an angle which is not orthogonal to the plane of the wake. 
Then, relative to photons emitted at the central point on 
the photon path, photons emitted from the highest point 
and the lowest point obtain a relative Doppler shift of 



61^ „ . ^Hw 

— = 2sm(6')tane' 

c 



(20) 



where H is the expansion rate of space and w is the wake 
width computed in the previous section, both evaluated 
at the redshift z when the photons are emitted. The 
angle 9 is indicated in Figure 2. It is the angle relative 
to the vertical to the wake. As a consequence of the 
normalization of 0(i^) we hence find 



4i{iy) = — for ve[iyio 

dv 2 



vio + — ] 



(21) 
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FIG. 2: Photons reaching us from a particular direction given 
by the angle 9 to the vertical to the wake are emitted at differ- 
ent points in the wake. Three such points are indicated in the 
Figure - the central point and two points on the bottom and 
top of the wake, respectively. Since the wake is undergoing 
Hubble expansion in its planar directions, there is a relative 
Doppler frequency shift Si/ between the 21cm photons from 
the top and the center of the wake. 



and — otherwise. 

In addition to its role in determining the brightness 
temperature, the frequency shift Sv is important since 
it determines the width of the 21cm signal of strings in 
the redshift direction and is hence central to the issue 
of observability of the signal. From ([20|) and ([8]) (and 
setting the 2sin9tan0 factor to one) we find that 

ly 15 

~ 3x 10-5(GAi)6(«s7.), (22) 

where in the second line we have used the cosmological 
parameters mentioned at the end of the introductory sec- 
tion and inserted the representative redshifts Zi + 1 = 10^ 
and z{t) + l = 30. 

Taking the formula P7|) for the brightness tempera- 
ture, inserting the results (fTSi) for the optical depth and 
(PU]) and (pij) for the line profile, we get 



STbii^) = 2- 



1- 



1/2 



xn%{to) 



^wake 



(to) 



bg 



(to) 



x(2sin2(e'))-i(l + z)i/^ 



(23) 



where r2,„ is the fraction of the critical energy den- 
sity which is in the form of matter. Here the ratio of 
the density inside the wake to the background density, 
nj^'}'"^ /n'^j, is approximately 4. We have rescaled the 
Hubble constant and the background HI density to its 
current values Hq and 'n'Jfj(to), respectively, and made 

use of H{z) = Hoi}m^{l + zf"!"^ in the re-scaling. Note 
that the width of the wake has cancelled out between 
the HI column density and the line profile. The width of 



the wake, however, has not disappeared completely from 
the calculation since it determines the wake temperature 
Tjf, and since it yields the width of the signal in redshift 
direction. 

Taking the values A^q = 2.85 x lO^^^ s^^, = 
0.068 K, = 73 km s"^ Mpc"\ v^q = 1420 MHz, 



fib = 0.042, = 0.26, 2sin^6' = 1, eq. 



becomes 



8T^{v) = [0.07 K] 



1 



(24) 



The collision coefficient Xc is dominated by the 
hydrogen-hydrogen collisions because of the very low 
fraction of free electrons. It is given by [6l| 



^10 



(25) 



where k^q is the de-excitation cross section and is ap- 
proximately 2.7x 10~^^cm''s^^ at a wake gas temperature 
of Tk = 20K [6l!| . This temperature corresponds to the 
parameter values we took at the end of Section 2, i.e. 
(G^)g = 0.3, (vsjs)^ = 1/3, 2, = 10^ and z + 1 = 30. 
Remembering that the density n inside the wake is four 
times the background density ribg = 1.9 x 10^^cm^'^(l -|- 
z)^ we find for a redshift 1 + z = 30 that Xc — 0.16 and 
hence (from 6Ti,{v) ~ —160 mK. For a formation 



redshift Zi = z^q = 3200 corresponding to the time of 
equal matter and radiation the corresponding tempera- 
ture is 6Ti,{fj,) ~ —15 mK. Note that the dependence 
of the above result on the cosmic string tension n enters 
only through the wake temperature. Note also that in the 
large G/i limit (G/i ^ 10~^), the brightness temperature 
approaches a constant value of close to 400mK. 

There is a critical value of the string tension (which 
depends on the redshift) at which the cosmic string sig- 
nal changes from emission to absorption. This value is 
determined by Tx = Tj which yields 



(Gm)^ ^ 0.1{vsl. 



.-2 (^+1)' 



Z, + 1 



0.4 



Zi + 1 



(26) 



for the value (vs^s)"^ = 1/3 which we are using through- 
out, and this corresponds to {G^)q ~ 0.6 if we use 
{zi + 1) = 10'^ and insert our reference redshift of 
z + 1 = 30. At a redshift of z -I- 1 = 20 the critical value is 
{Gfj,)e ~ 0.4. For a formation time corresponding to the 
redshift of equal matter and radiation the critical value 
of G/i is a factor of two lower. 

For lower values of the tension the cosmic string signal 
in the 21cm maps shifts from an emission signal to an 
absorption signal. However, two effects must be taken 
into account when discussing the visibility of the string 
signal. The first is the fact that at very small values 
of G/i our calculation of the wake temperature breaks 
down. This occurs when the resulting wake temperature 
computed in ()13|) is smaller than the gas temperature 
obtained by taking the background gas temperature and 
computing its increase under adiabatic compression to 
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the overdensity of the wake. In this case it is no longer 
justified to take the initial conditions ([5]). At redshifts z 
below 150 the cosmic gas is cooling adiabatically as 

Tg = 0.02K(1 + zf (27) 

because at this point Compton heating of the gas by the 
CMB is negligible Setting Tk = 2.5 x Tg (the factor 
2.5 being due to the fact that for adiabatic compression of 
a mono-atomic gas by a factor of 4 in volume, one expects 
a temperature increase by a factor of 4^/'^ '-^ 2.5) yields 
(for our standard values of 1 + = 10"^ and [vs^sY = 
1/3) 

(G/i)6 =i 3 X 10-^(1 + z)3/2 (28) 

which equals {Gfi)e — 0.5 at a redshift of 1 + z = 30. 
For a formation redshift Zi — z^q the factor 3 in (|28p is 
replaced by 2. 

For values of Gfi smaller than the one given in (pS)) , the 
effects of thermal pressure start to become important, 
and this will effect both the density and temperature 
structure in the wake. 

Our results are illustrated in Figure 3 which shows a 
comparison of the temperatures relevant to the above 
discussion and their dependence on redshift. The ver- 
tical axis is the temperature axis, the horizontal axis is 
inverse redshift. The magenta (dashed) line shows the 
CMB temperature T^, the orange (dotted) line is the gas 
temperature in a wake 2.5Tg after adiabatic contraction. 
The two solid lines with the positive slope represent Tk 
for two different values of Gfi. The upper curve is for 
{Gfi)e — 1, the lower curve for {Gfi)e = 0.3. Along a 
fixed G/Lt curve, the 21cm signal of a cosmic string is an 
emission signal above the Tj curve and an absorption 
signal below it. Above the 2.5Tg curve, the tempera- 
ture of the gas inside the wake is well approximated by 
the equation (jl3p . below it the initial thermal gas tem- 
perature effects dominate and the wake temperature will 
follow the 2.5Tg curve. The wake temperature curves are 
for a formation redshift of -I- 1 = 10'^. 

The second issue is that the string signal must be com- 
pared to what would be seen if the wake region were 
replaced by unclustered neutral gas. In this case, the 
absorption signal has a temperature given by ([M]) with 
Tk replaced by Tg and the collision coefficient computed 
with Tg instead of Tk- Due to the overdensity of the 
wake, the signal of a string wake should persist. The 
overdensity in the string wake by a factor of 4 will lead 
to an enhancement in the magnitude of the brightness 
temperature by a factor of 16 compared to what would 
be seen if the wake region were replaced by unclustered 
gas. One factor of 4 comes from the fact that the collision 
coefficient is proportional to the gas density, the second 
from the fact that the optical depth is also proportional 
to the density. In addition, the temperature ratios in ([M)) 
are different. 

However, for values of G/i so small that Tk ^ Tg there 
will no longer be any shock and hence no well-defined 
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100 60 50 30 20 10 

z 

FIG. 3: The redshift dependence of the relevant tempera- 
tures. The majenta dashed line is the CMB temperature Tj, 
the orange dotted line is 2.5Tg, the gas temperature after 
adiabatic contraction by a factor 4 in volume, and the green 
(dot-dashed) line is the Tg curve. The two solid curves with 
positive slope represent Tk for two different values of Gfi, 
{Gfi)6 = 1 in the case of the upper curve and {Gfj.)^ = 0.3 
in the case of the lower curve (in both cases for a formation 
redshift Zi + 1 = 10^). Following one of the solid curves, we 
see that the 21cm string signal is an emission signal above the 
majenta curve and an absorption signal below it. The wake 
temperature is given by Tk only above the orange curve. For 
higher redshifts, the initial gas temperature dominates the fi- 
nal temperature of the gas inside the wake, and the latter fol- 
lows the orange curve. Once Tk drops below the green curve, 
there will no longer be any shock and our efi'ect disappears. 



wake region with overdensity 4. This will occur at a value 
of G^ which is smaller than the limit in (|28p by a factor 
of ^2.5. Note also that the thickness in redshift space 
(discussed below) of the signal decreases as G/i decreases, 
and hence improved sensitivity will be required to detect 
the signal. 

The bottom line of the above analysis is that at a red- 
shift of z + 1 = 30, strings with tensions below the current 
observational limit are predicted to be visible in absorp- 
tion in 21cm surveys. Strings with larger tensions would 
yield an emission signal. The critical value ((26)) of G/i at 
which the emission signal turns into an absorption signal 
decreases as z decreases, so that strings with tensions at 
the current observational limit would become visible in 
emission at a redshift of below 20. However, the value 
of z cannot be smaller than that corresponding to reion- 
ization since we have not considered how UV scattering 
affects the spin temperature. 

The 21cm signature of a cosmic string wake has a dis- 
tinctive shape in redshift space. The string signal will be 
a wedge-like region of extra 21cm emission. The signal 
will be wide in the two angular directions, and narrow in 
redshift direction ^. The projection of this region onto 



The signal at a fixed redshift would be extended only in one 
angular direction. 
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the plane corresponding to the two angles in the sky (we 
are working in the limit of small angles and thus can 
use the flat sky approximation) corresponds to the pro- 
jection of the wake onto the observer's past light cone. 
The orientation of the wedge in redshift space is given 
by the orientation of the wake relative to the observer's 
light cone. This is illustrated in Figure 4. The left panel 
is a space-time sketch of the geometry (with two spatial 
directions suppressed). The wake is created by a string 
segment which starts out at time ti at the position xi and 
which at time 2ti (roughly one Hubble time step later) 
has moved to the position X2 (the arrow on the segment 
connecting the two events gives the direction of the ve- 
locity of the string segment. For the orientation of the 
string chosen, the past light cone intersects the "back" 
of the string wake at a later time than the front. Hence, 
the 21cm signal has a larger red-shift for photons from 
the tip of the wake than from the back side of the wake. 
The width of the wake vanishes at the tip and increases 
towards the back. Hence, in redshift space the region of 
extra 21cm emission has the orientation sketched on the 
right panel of Figure 4. The planar dimensions of the re- 
gion of extra emission are in the angular directions, and 
their sizes are the angles corresponding to the comoving 
area (see ^) 



ClU{zi + 1) X tiVs^siZi + 1) 



(29) 



The amplitude of the emission signal depends on the ori- 
entation of the wake relative to us (through the redshift 
when our past light cone intersects the various parts of 
the wake). This is not correlated with the direction of 
string motion. The width in redshift space, on the other 
hand, depends on the emission point on the wake - the 
width is larger in the back and approaches zero at the 
tip, as given by Si/ (see [20)) . 



IV. CONCLUSIONS AND DISCUSSION 

In this paper we have calculated the 21cm signal of a 
single cosmic string wake. We found that wakes leave 
a characteristic signal in 21cm surveys - wedge shaped 
regions of extra emission or absorption, depending on 
the value of the string tension. There is an emission 
signal as long as the string tension exceeds the critical 
value given by (l26l) and as the tension increases the tem- 
perature approaches the value given by 200mK, a value 
larger by almost two orders of magnitude compared to 
the backgrounds from standard cosmology sources which 
the Square Kilometer Array (SKA) telescope project is 
designed to be sensitive to [63]. For values of below 
the critical value given by ([25)) . the cosmic string wake 
signal changes to an absorption signal. This will persist 
down to a value of G/i a factor of y/2.5 smaller than the 
value given in (|28p. at which point the material which is 
gravitationally attracted to the wake will no longer un- 
dergo a shock and hence our analysis ceases to be appli- 
cable. The width of the wedge in redshift space depends 
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FIG. 4: Geometry of the 21cm signal of a cosmic string wake. 
The left panel is a sketch of the geometry of the wake in 
space-time - vertical axis denoting time, the horizontal axis 
one direction of space. The string segment producing the 
wake is born at time ti and travels in the direction of the ar- 
row, ending at the position X2 at the time 2ti. The past hght 
cone of the observer at the time to intersects the tip of the 
wake at the time S2, the back of the wake at time si. These 
times are in general different. Hence, the 21cm radiation from 
different parts of the wake is observed at different red-shifts. 
The resulting angle-redshift signal of the string wake shown 
in the left panel is illustrated in the right panel, where the 
horizontal axis is the same spatial coordinate as in the left 
panel, but the vertical axis is the redshift of the 21cm radi- 
ation signal. The wedge in 21cm has vanishing thickness at 
the tip of the wedge, and thickness given by Si/ at the back 
side. 



on the string tension via (|20[) . The width increases from 
zero at the tip of the string to the value given by ([20)) at 
the midpoint of the wake. 

The sensitivity of 21cm surveys to cosmic strings is 
best at the lowest redshifts sufhciently higher than the 
redshift of reionization. This is due to the fact that the 
string wakes keep increasing in width as a function of 
time, leading to an increasing temperature of the gas in 
the wake and hence to a higher 21cm signal. At the same 
time, the wakes become more stable towards thermal dis- 
ruption. 

As mentioned above, the brightness temperature shift 
which we predict will be large compared to that from the 
surrounding intergalactic medium. The fractional fre- 
quency width of the signal is given by ([22)) and is near 
the limit of the frequency resolution of radio telescopes 
for values of G/x close to the current upper bound. For 
example, the Square Kilometer Array will have a frac- 
tional spectral resolution of the order [U of 10^*. In 
both angular directions, the signal will cover a rectangle 
whose side length is given by the comoving Hubble radius 
at the redshift Zi, which for redshifts close to recombina- 
tion corresponds to a scale of about one degree. Thus, 
the signal of cosmic strings which we predict is clearly 
a possible target for future radio telescopes. For exam- 
ple, the angular resolution of the SKA is designed to be 
smaller than 0.1 arcsec, and the frequency range of the 
SKA will allow the detection of the 21cm signal up to a 
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redshift of 20. 

We have focused on the signal of a single wake laid 
down at a time t > t^q. The reader may worry that 
the signal of such a wake is masked by the "noise" due 
to wakes laid down at earlier times. Even though the 
baryons will stream out of these wakes due to their cou- 
pling with the photons, the dark matter wakes persist. 
For t < teq the dark matter wakes do not grow in thick- 
ness until t = teq. After that, the thickness will begin 
to grow, and at t > tree the baryons will start to fall 
into the dark matter potential wells. However, the width 
of such wakes (averaged over the length of the wake) is 
smaller than that of wakes at t = teq. In addition, due to 
their large number, on length scales of wakes laid down at 
t > teq, the earlier wakes will act like Gaussian noise. The 
coherent signature in position space of such a "late" wake 
can be picked out of the Gaussian noise even if the am- 
plitude of the Gaussian noise is comparable to or slightly 
larger than the amplitude of the signal (as measured in 
terms of the contribution to the power spectrum). This 
has been studied in the context of picking out the signa- 
tures of late strings in CMB temperature maps in [43 - 1441 , 
and we expect similar conclusions to hold here. 

Let us end with a brief comparison of our work with 



that of [58] and [539] who also considered 21cm signals of 
cosmic strings. What sets our work apart is that we focus 
on the specific position space signature of wakes rather 
than on the (Fourier space) power spectrum. In comput- 
ing a power spectrum, one loses the information about 
the non-Gaussianities in the distribution of strings. It 
is these non-Gaussianities which most clearly distinguish 
the predictions of a string model from a model with only 
Gaussian fluctuations. Therefore, the sensitivity to the 
presence of cosmic strings will be much higher in a study 
like ours compared to what can be achieved by only com- 
puting power spectra. 
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